Customer No. 20350 *' 

TOWNSEND and TOWNSEND and CREW LLP « | 
Two Embarcadero Center, 8 th Floor ^ [ 

San Francisco, California 941 1 1-3834 I 
(415)576-0200 JJj 



Attorney Docket No. 



17003-23. 00US 



"Express Mail" Label No. . 
Date of Deposit: 



EL265823664US 



I hereby certify that this is being deposited with the United States 
Postal Service "Express Mail Post Office to Addressee" service 
under 37 CFR 1 .10 on the date indicated above, addressed to: 



Assistant Commissioner for Patents 
Washington, D.C, 20231 



ASSISTANT COMMISSIONER FOR PATENT&^= 
BOX PATENT APPLICATION 
Washington, D.C. 20231 

Sir: 

Transmitted herewith for filing under 37 CFR 1.53(b) is the 
[ X ] patent application of 
[ ] continuation patent application of 
[ ] divisional patent application of 
[ ] continuation-in-part patent application of 

Inventor(s)/Applicant Identifier: Fredric S. Young 



For: METHOD AND APPARATUS FOR DESCRIBING AND SIMULATING COMPLEX SYSTEMS 
[ ] This application claims priority from each of the following Application Nos./filing dates: 



September 23> 1999 




[ ] 



ifiiclosed 

*] 
*] 

-X] 

on 



II] 



the disclosure(s) of which is (are) incorporated by reference. 
Please amend this application by adding the following before the first sentence: "This application is a 

continuation-in-part of and claims the benefit of U.S. Application No. 60/ , filed 

which is incorporated by reference." 

are: 

_U_page(s) of specification 
2 page(s) of claims 
1 page of Abstract 

8 sheet(s) of [ ] formal [X] informal drawing(s). 
An assignment of the invention to Chroma Graphics. Inc. 



] continuation [ ] 
. the disclosure of 



A signed Declaration & Power of Attorney 

A verified statement to establish small entity status under 37 CFR 1.9 and 37 CFR 1.27 is enclosed 
Information Disclosure Statement under 37 CFR 1.97 
Form PTO-1449 
One reference 



U FOR: 


NO. FILED 


NO. EXTRA 




TOTAL 


10 -20 


*o 


CLAIMS 






INDEP. 


2 -3 


*o 


CLAIMS 






f 1 MULTIPLE DEPENDENT CLAIM PRESENTED 



SMALL ENTITY 



OTHER THAN 
SMALL ENTITY 



* If the difference in Col. 1 is less than 0, enter "0" in 
Col. 2. 

Please charge Deposit Account No. 20-1430 as follows: 
[ X ] Filing fee 



RATE 


FEE 


OR 


RATE 


FEE 




$380.00 


OR 




$760.00 


x $9.00 = 


$0.00 


OR 


x $18.00 = 




x S39.00 = 


$0.00 


OR 


x $78.00 = 




+ $130.00 = 




OR 


+ $260.00 = 




TOTAL 


$380.00 


OR 


TOTAL 





$380.00 



[ X ] Any additional fees associated with this paper or during the pendency of this application. 
[ ] The issue fee set in 37 CFR 1 . 1 8 at or before mailing of the Notice of Allowance, pursuant to 37 CFR 1.311(b) 



L 1 



A check for $ 



is enclosed. 



2 extra copies of this sheet are enclosed. 



Respectfully submitted, 

TOWNSEND and TOWNSEND and CREW LLP 



Telephone: 
(415) 576-0200 



Facsimile: 
(415) 576-0300 



Kenneth R. Allen 
Reg No.: 27,301 
Attorneys for Applicant 



PA 3023340 v1 



Attorney Docket No.: 17003-23.00US 



VERIFIED STATEMENT (DECLARATION) CLAIMING SMALL ENTITY STATUS 
(37 CFR 1.9(f) & 1.27(c)) - SMALL BUSINESS CONCERN 



Applicant or Patentee: Fredric S. Young 

Application or Patent No.: 



Filed or Issued: Herewith 



Title: METHOD AND APPARATUS FOR DESCRIBING AND SIMULATING COMPLEX SYSTEMS 



I hereby declare that I am: 



[ ] the owner of the small business concern identified below: 

[x] an official of the small business concern empowered to act on behalf of the concern identified below. 



Name of Small Business Concern: Chroma Graphics, Inc. 



Address of Small Business Concern: 577 Airport Blvd., Suite 730 



Burlingame, CA 94010 



I hereby declare that the above-identified small business concern qualifies as a small business concern as defined in 13 CFR 121.12, and reproduced in 37 CFR 
1.9(d), for purposes of paying reduced fees to the United States Patent and Trademark Office, in that the number of employees of the concern, including those of 
its affiliates, does not exceed 500 persons. For purposes of this statement, (1) the number of employees of the business concern is the average over the previous 
fiscal year of the concern of the persons employed on a full-time, part-time or temporary basis during each of the pay periods of the fiscal year, and (2) concerns 
are affiliates of each other when either, directly or indirectly, one concern controls or has the power to control the other, or a third party or parties controls or has 
the power to control both. 



I hereby declare that rights under contract or law have been conveyed to and remain with the small business concern identified above with regard to the invention, 
entitled METHOD AND APPARATUS FOR DESCRIBING AND SIMULATING COMPLEX SYSTEMS by inventor(s) Fredric S. Young described in: 



«? [x] the specification filed herewith; 

13 [ ] Application No. , filed ; 

jp' [ ] Patent No. , issued . 

^If the rights held by the above identified small business concern are not exclusive, each individual, concern or organization having rights in the invention is listed 
: ^below* and no rights to the invention are held by any person, other than the inventor, who would not qualify as an independent inventor under 37 CFR 1.9(c) if 
. inthat person made the invention, or by any concern that would not qualify as a small business concern under 37 CFR 1.9(d), or a nonprofit organization under 37 
^CFR 1.9(e). 

*NOTE: Separate verified statements are required from each named person, concern or organization having rights to the invention averring to their 
v status as small entities. (37 CFR 1.27) 

:^,Name: 
Addres s : 

y [ ] Individual [ ] Small Business Concern [ ] Nonprofit Organization 

Zl acknowledge the duty to file, in this application or patent, notification of any change in status resulting in loss of entitlement to small entity status prior to paying, 
yor at the time of paying, the earliest of the issue fee or any maintenance fee due after the date on which status as a small entity is no longer appropriate. (37 CFR 
1 .28(b)) 



I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are believed to be true; and 
further that these statements were made with the knowledge that willful false statements and the like so made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States Code, and that such willful false statements may jeopardize the validity of the application, any patent issuing 
thereon, or any patent to which this verified statement is directed. 



Name of Person Signing: Arch F. Meredith HI 



Title of Person Other than Owner: Executive Vice President 



Address of Per/on Eigning y''T^r- 577 Airport Blvd., Suite 730 




Burlingame, CA 94010 

Signature V^Mj^^r ' ^Ou. ">> Date 




PA 3023066 v1 



Page 1 of 1 



Attorney Docket No.: 17003-23.00 



PATENT APPLICATION 

METHOD AND APPARATUS FOR DESCRIBING AND 
SIMULATING COMPLEX SYSTEMS 



Inventor: 

Fredric S. Young 
128 Lyell Street 
Los Altos, C A 94022 
a citizen of the USA 



Assignee: 

Chroma Graphics, Inc. 
577 Airport Blvd., Suite 730 
Burlingarne, CA 94010 



Entity: Small 



TOWNSEND and TOWNSEND and CREW LLP 
Two Embarcadero Center, 8 th Floor 
San Francisco, California 941 1 1-3834 
(650) 326-2400 



1 

PATENT 

Attorney Docket No.: 17003-23.00 



METHOD AND APPARATUS FOR DESCRIBING AND 
SIMULATING COMPLEX SYSTEMS 

5 

BACKGROUND OF THE INVENTION 
This invention relates to symbolic representation of data for both static and 
dynamic analysis and manipulation. In particular, the invention relates to methods and 
apparatus for modeling n-dimensional complex systems and networks, including 

10 biological systems, social systems, geological formations and processes and simulations 
thereof This invention employs a visual calculus of complex systems alluded to in the 
Ph.D. dissertation of the present inventor entitled "Determination and Stabilization of the 
Bacterial Growth Rate " by Fredric S. Young, University of Michigan, 1977. 

Prior work including the aforementioned dissertation of the present 

1 5 inventor failed to suggest anything beyond a simple modeling of simple cellular processes 
in simple single prokaryotic cell types which include bacteria, and prior work has failed 
to describe or suggest how cells might interact in multicellular organisms. The 
dissertation of the present inventor was a theory of the computational processes of 
bacterial cells in natural and synthetic environments. This dissertation was an early effort 

20 in what has now come to be known as bioinformatics. Since all multicellular organisms 
important in medicine and physiology contain the much more complicated eukaryotic 
type cells, the prior work was not applicable to models other than simple bacterial cells. 

In a parallel development, Per Bak at Brookhaven National Laboratories 
proposed a general model for complex systems to explain the ubiquitous occurrence of 

25 fractal structures and fractal (1/f) noise in a wide variety of physical and other natural 
systems. 

It has been observed that certain non -equilibrium processes cannot be 
described and analyzed with sufficient mathematical clarity with current mathematical 
tools. Efforts have been made in recent years to develop the mathematics of nonlinear 
30 systems using nonlinear dynamics and complexity theory. An interesting and major 

lesson learned from the dissertation research of the present inventor and the later research 
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in non-linear dynamics is that there are simple alternatives to conventional differential 
equation based simulation that can capture the essence of a complex system in a greatly 
simplified or "toy" model. 

What is needed are techniques and devices to exploit these discoveries for 
description and ultimately simulation in some of the most economically significant 
applications and problems in geology, biology and economics. Simulation models can 
then be substituted for laboratory and field research to guide diagnosis and therapy 
development in medicine, data processing and decisionmaking affecting the acquisition 
and development of natural resources. 



SUMMARY OF THE INVENTION 
According to the invention, a method for description and simulation based 
on organizing data into maps of invariants, the invariants being points of entropy balance 
in a system of interest which is either in a stationary state or in a transitory disturbed state. 

1 5 The method includes identifying invariants in the system of interest by identifying 
primary sources and sinks of energy, identifying secondary energy sources and sinks 
coupled to the primary sources and sinks, and coupling all such sources and sinks into a 
network of transformations organized around nodes of those sources and sinks 
corresponding to the invariants, each of the nodes being characterized by a locally defined 

20 principle of balanced self-organization in a system with both a conservation law and 
energy dissipation. Such a system becomes "organized" upon achievement of a critical 
rate of entropy flux into the environment. Associated with each invariant are response 
rates related to energy transfer rates into and out of the invariants. 

The invention is based on the discovery that all systems subject to input of 

25 any source of energy are rendered stable where the conserved quantity is at a local "angle 
of repose," that is, where all input rates and output rates are balanced with respect to 
energy input and dissipation. Systems organized in this manner exhibit either first or 
second non-equilibrium phase transitions. Systems with second order phase transitions 
have critical points whose properties are mathematically related to the critical points 

30 found in select rare equilibrium systems that undergo second order phase transitions, such 
as systems that have a point separating three phases, as shown in a physical phase 
diagram. (Examples include pressure/temperature effects on water and on carbon 
dioxide.) Invariants are conserved ratios reflecting the local angles of repose that result 
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from the conservation of quantities in a dissipative system. These invariants associated 
with critical ratios can be described using the mathematics of percolation to describe the 
dynamics at the critical point. 

Experimental studies of self-organization have shown, in contrast to the 
5 suggestion of Bak, that generic non-equilibrium self-organization is most likely to 

organize systems at first order and not at second order phase transitions. Self-organization 
at second order phase transitions can be achieved by incorporating extremum principles 
allowing selection for maxima and minima in energy dissipation and optimization. A 
system organized at a first order phase transition can be brought to the critical point of a 

10 second order phase transition by including feedback for optimization. The subset of 
systems containing second order critical points includes both living and non-living 
systems. According to the invention, the critical points in living systems are stabilized by 
a unique form of double reciprocal feedback. 

According to one application of the invention, a method is provided for 

1 5 controlling the engineering of a complex system by creating a database based on a series 
of steps beginning with a description of a structure, organism or system to be created 
which was obtained by tracking flow of energy and transformation of "atomic" 
(undivisible) elements into complex structures along a reaction chain; then, given such a 
description of a structure, recognizing that it is a self-organizing model of the complex 

20 system which could have either a first order or a second order phase transition and that 
the scope of behaviors of the system is limited by the description. (A first order phase 
transition has no critical point, whereas a second order phase transition has a critical point 
in which at least three phases co-exist: i.e., solid, liquid and gas.) This is responsible for 
the phenomenon of universality wherein a great number of systems exhibit common 

25 behavior as evidenced in power laws describing energy scaling in systems that are valid 
over a huge range of scales. Knowing that non-equilibrium systems can only become 
balanced in a system with a phase transition of the first order or of the second order, and 
knowing that all living systems are characterized by biological regulatory mechanisms 
that are stabilized at the second order phase transition point using a double feedback 

30 mechanism, it follows that all living systems have a pattern of allowable behaviors that 
can be predicted and thus designed. 

The present invention is a nontrivial extension of the Ph.D. work of the 
present inventor to n-dimensional analysis and from the growth of bacteria to homeostasis 



in all biological systems as well as stationary states in all systems not in equilibrium (that 
is, systems subject to input and dissipation of energy). 

In a specific embodiment as an example of a workflow description, each 
invariant may be modeled as a storage element having an activated or energized state and 
an inactive or unenergized state with an internal mechanism for transitioning over a time 
scale between the two states, with unidirectional flow of energy through every storage 
element. The invariant is mechanically modeled as two variable-rate one-way valves 
connected through an accumulator, wherein the net flow of all valves in a system must be 
zero in the stationary (undisturbed) state. This condition of net flow of N valve sets 
equaling zero can thus be satisfied only when each valve set in the network is at a local 
critical point. In other words, at each invariant there is a fixed relationship between the 
flow of each valve set in the network and the percentage of the energy stored in the 
energy storer. This relationship can be completely defined whenever a dissipative system 
can be associated with a conserved quantity. 

Organizing a simulation around invariants produces a set of constraints 
which allows the development of minimally complex representations of any system of 
interest. 

One of the advantages of the invention is the provisions a nearly absolute 
reference framework for organizing any additional data into the simulations. In the 
examples of homeostatic regulation in biology, the inventive method of organization of a 
simulation provides a method of vastly streamlining the work involved in the human 
genome project. 

The invention will be better understood by reference to the following 
detailed description in connection with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a multiscale spatial static dataset as might be derived by applying 
pattern recognition techniques to geological data and illustrates three levels of scale in the 
indexing of the tree and three frequencies of branching events at each of the levels. 

Fig. 2A-Fig. 2F is an illustration of the first five levels of a fractal of the 
type known as a space filling curve illustrating the case where levels are homogeneously 
distributed. 



Fig. 3 is an illustration of the first four levels of a fractal of the type known 
as a space filling curve illustrating the case where levels are not homogeneously 
distributed. 

Fig. 4 is multilevel fractal curve drawn as a one dimensional temporal 

dataset 

Fig. 5 A and 5B is an illustration of a result of a process according to the 

invention. 

Fig. 6 is an illustration of a tool employed in a first step in a process 
according to the invention. 

Fig . 7 is a flow chart of the method for fractal modeling according to the 

invention. 

Fig. 8 A and 8B is a flow chart of workflow process development 
according to the invention. 

Fig. 8 is an apparatus according to the invention. 

Fig. 9 is a multidimensional network of superposed multilevel objects 
according to the invention. 

Fig. 10 is a diagram of a multidimensional network of a plurality of nodes 
of the type of Fig. 9. 

DESCRIPTION OF SPECIFIC EMBODIMENTS 
In order to understand the invention, it is useful to define the underlying 
elements. Referring to Fig. 1, there is illustrated a graphical representation of a typical 
multiscale spatial static dataset 10 as might be derived by applying pattern recognition 
techniques to geological data. It illustrates three scale levels 12, 14, 16 in the indexing of 
the tree. Each scale level has three branching event frequencies. ( At each increase in 
scale, to each source line 18, 19 is added a corresponding triangle 20, 21 resulting from 
events at the third harmonic of the source line.) This dataset 10 is a multiscale fractal, 
meaning that the fractalization process is not applied homogeneously at each scale. For 
the purposes of explanation, the pattern shown is totally regular, each level is of the same 
exponential. An example of an irregular pattern would be the boundary of a seacoast 
viewed at different scales. However, each level would be covered by a range of 
exponentials. The dataset 10 has the same characteristics of any dataset derived from an 
appropriate pattern recognition system which can yield labeled textures. An example of a 
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suitable pattern recognition system is described in U.S. Pat. Application Serial Number 
09/070,1 10 filed 4/29/98. Other pattern recognition systems may provide similar results. 

Fig. 2A-Fig. 2F are illustrations of the first five levels 22, 24, 26, 28, 30 of 
a fractal of the type known as a space filling curve illustrating the case where levels are 
5 homogeneously distributed. This illustration of a multilevel space-filling curve. Multi- 
level space-filling curves can be combined with other types of curves as hereinafter 
explained. The process of progressing through each of the levels is called fractalization. 

Fig. 3 is an illustration of the first four levels 32, 34, 36, 38 of a fractal 40 
of the type known as a space filling curve illustrating the case where levels are not 
10 homogeneously distributed. The process of producing this fractalization is a variant of the 
type shown in Fig. 1, operating on each line segment of Fig. 2A-2F. The output of a 
pattern recognition system yielding labeled textures would be as illustrated in Fig. 3. 

Fig. 4 illustrates a multilevel fractal curve 42 drawn as a one-dimensional 
temporal dataset wherein the local density of the curve is a measure of the level of fractal 
15 branching of the type of process shown in Fig. 2A-2F. Thus the fractal 42 is a 

combination of a one dimensional multilevel temporal dataset and a space filling curve. 

Figs. 5 A and 5B are an illustration of a result of an environmental process 
44 which is a superposition of temporal processes 42 and the datasets previously 
described. This illustrates how objects defined as labeled textures may process 
20 throughputs, and it represents a graphical description of a physical process analyzed 

according to the invention. The multiscale nature of the process should be evident from 
the illustration, which shows the magnification in Fig 5B of one of the processes 47 as 
having the same characteristics and structure as the process 42 in Fig. 5A from which it is 
magnified. 

25 Fig. 6 is an illustration of a system 46 for producing the description or 

simulation 45 according to the exact transformation process 48 of the invention. The 
inputs to the process 48 are the three dimensional fractal dataset DS1 10 of Fig. 1 and the 
one-dimensional temporal dataset DS2(t) 42 of Fig. 4. A computer program operative 
according to the inventive process 48 would produce the four dimensional simulation set 

30 DS3 or description 45 in accordance with the invention. Not shown but which should be 
understood is that the output product 45 can produce a second completely orthogonal one- 
dimensional temporal dataset which may be the basis of feedback to each dimension of 
the input datasets 10 and 42, changing both the temporal and organizational 



characteristics of the transformation process 48. (This feedback is not illustrated in Figure 
6, for simplicity, because it would be in a dimension not readily illustrated in the 
dimensional depiction of Fig. 6. However, such depiction could be modeled in grayscale 
or color.) This second one-dimensional dataset can be a mapping of other higher 
dimensional information. A pair of related four dimensional datasets of space and time 1 
and space and time 2 would constitute a six dimensional dataset. 

Fig . 7 is a flow chart of the process 48 for fractal modeling according to 
the invention. The process begins by a gross form of data compression, namely, the 
segmenting the data into different textures (Step A). Textures are patterns or 
combinations of patterns. Examples are identifiable statistically repeated patterns as 
might be found in rock formations, chains of DNA and the like. 

The textures are then labeled (Step B) for appropriate identification. The 
labels serve as a tags for the compressed data resulting from preliminary analysis based 
on pattern/texture recognition techniques. 

The process is furthered by defining the system and its environment, as 
well as the boundaries between the system and the environment (Step C). The system is 
the texture to be considered, and the environment is everything that affects the system 
which is necessary to make the model a closed model where the system is an open 
system. 

Thereafter comes a description of the workflow of the system (Step D). 
This workflow must be in terms of sources and sinks of energy and of raw materials, as 
hereinafter noted in Figs. 8 A and 8B. 

Referring to Figs. 8A and 8B, first, a list of the energy and elementary 
materials of the system is developed (Step E). 

Then, for each member of the list, the points of entry into the system from 
the environment are identified (Step F). The points of entry are subsets of the labeled 
textures. 

Then, for each point of entry, roots of the points of entry are traced 
through the system, and points of further transformation, if any, are identified (Step G). 
Points of further transformation can be fractal level changes, changes in scale, changes in 
texture, changes in content as by splitting or joining or rearrangements and reordering of 
content. 
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When that is done (Step H), the nodes are characterized by their inputs, 
outputs and transformations (Step I). Each transformation is described as a process in a 
workflow diagram. 

When that is done (Step J), the relative rate of each process is catalogued 
as either balanced, unbalanced fast or unbalanced slow, in a three state system (Step K). 
For each level there is a unique temporal behavior. When there is a stationary state 
defined over the whole system by an interrelated state of balance, the stationary state is 
defined relative to a lowest interval of time that can be analyzed. If there is an attempted 
analysis with a finer increment of time, for example, fluctuations away from the 
stationary state will be observed. At a second order phase transition, these fluctuations 
have no fixed scale. Since the set points are dynamic, the behavior in approach to repose 
from fluctuations is also without scale. This is analogous to the observation that 
avalanches can occur in all scales of piles of rubble, from sand to boulders. 

The stationary state is thus characterized statistically as zero error in time, 
i.e., it will average to zero error on a selected time scale. However, for the 
characterization in time there will also be a multilevel behavior of the pattern with 
deviations away from the stationary state. 

After the relative rates of the processes have been catalogued (Step L), 
each process is catalogued by level and frequency (Step M). In other words, a 
determination is made as to where each process is in the system and how many 
occurrences of the processes exist within the system. 

When that is done (Step N), the model is completed by mapping each 
process to a level and frequency with appropriate description in the level/frequency 
diagram (Step O). An example of a level/frequency diagram has been shown in Fig. 5. 
For each unique processor i (defined by a unique level and frequency), there is a 
workflow description of the process which is going on. This workflow description is 
embedded in the processor i, which is an object in an object oriented system. 

A key step according to the invention is the cataloguing of the relative 
rates (Step J). This involves tracking processes through various dimensions, levels and 
determining which side a process is relative to a locally-defined critical point. This 
important modeling step requires some external empirically-based input to effect. The 
selected choice of such a constant determines for example whether a system is self- 
perpetuating, increasingly oscillatory or decaying. Biological homeostasis is a 



particularly apt example of a self-perpetuating system where the system must be set at the 
self-perpetuating critical point. 

Fig. 9 is a diagram of a generic node representing the basic process 48 
according to the invention, and it best illustrates the basic process in all of its temporal 
manifestations. It is understood that this process is multidimensional and fractal in 
nature, according to the invention. Input and feedback can be from any fractal level or 
any associated dimension. The process comprises two stable coexistent states T x 50 and 
T 0 52 with first generalized activation 54 from T x 50 to T 0 52 and second generalized 
activation 56 from T 0 52 to Ti 50, and one or a plurality of energy and material inputs 58, 
60, 62 and an energy and material output 64 (which can always be represented as a single 
output). The generalized activations 54 and 56 are each a flow of energy and materials. 
The ratio of the (energy and material) population of state T x 50 to state T 0 52 is constant 
when the system is stable. For a living system, the process 48 represents homeostasis, that 
is, life in balance, or steady state. For a generalized physical system, the process 
represents stationary state. The process 48 is subject to internal self regulation, as 
hereinafter described, about a preselected local critical point, which serves to establish the 
local characteristics of the process and thus the ratio of stable states T]/T 0 which is 
defined as the steady state. When each local process is at steady state, a global system 
comprising the totality of local processes is also at steady state. That steady state is, 
according to the invention, defined as the self-organized critical state, and it is a critical 
point as found in a second order phase transition. If the global system cannot achieve 
criticality, then it is not a self-organized critical system. It may well self organize around 
a first order phase transition, but it is not a self organized critical system, since it lacks a 
critical point. 

The control system which establishes steady state is a feedback system 
with an error detector 68 for detecting deviations between the preselected ratio balance 70 
and the measured ratio balance 72. The error detector 68 controls an amplifier, or error 
control response subsystem 76, which in turn regulates inputs via valves 78, 80, 82 on 
inputs A, B and C, and, via a sequence time storage unit 84, a valve 86. (The sequence 
time storage unit 86 provides the time delay to assure that inputs and output are 
synchronized.) The preselected ratio balance 70 is established by an external critical 
point setter 74. The critical point is derived from the boundary conditions, structure and 
entropy considerations. Given a set of boundary conditions between the internal 
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mechanisms and the environment, the critical point can be calculated using graph theory 
and thermodynamics, i.e., the Second Law of Thermodynamics. Alternatively, the 
critical point can be empirically derived by comparing a real system and its corresponding 
simulation model. 

Fig. 10 is a diagram of a multidimensional network 49 of a plurality of 
nodes (processes) 48 in one dimension and processes 148 in another dimension of the 
type of Fig. 9. Feedback paths 51, 53 may be across dimensions. 

As an example of a model of a system according to the invention, consider 
the complexity of eukaiyotic cell growth. It is postulated from the model according to the 
invention that there must be an additional level of regulation to balance the growth 
allowed by the nutritional resources with the organism's need for those cells in the overall 
system. Cancer is an example of uncontrolled growth of particular classes of cells. When 
the first cancer-related gene was identified, it was shown to be in a protein that monitored 
the presence of growth factors and coupled this sensing to the allowable rate of cell 
growth. When this cancer gene was analyzed to determine if it was related to any known 
genes in bacteria, it was found to be related to the elongation factor. This relation was in 
fact found to be due to the same binding of GDP and GTP by the two proteins, namely 
the elongation factor and the cancer factor RAS. This can be modeled adequately by the 
present invention. In fact it can be used to model every known regulatory mechanism 
known to be relevant to cancer. 

As another example, cells are made of four classes of macromolecules, 
namely lipids, carbohydrates, nucleic acids and proteins. While these are important 
processes according to the definitional structure of the invention, they need not be 
modeled at a certain level and may well only be treated as input resources to a differently- 
indexed level of the system. The models of the cells can be subsumed into a subsystem 
such as an organ for purposes of macro simulation. 

It must be understood is that there is a thermodynamic critical point of 
phase transition associated each dimension of a system. Any system in a stationary or 
steady state has an inherent thermodynamic critical point related to a throughput of 
factors of the system. Within the limits of observation, if one can define a conserved 
quantity in the stationary or steady state of the system, then one can use a dynamic 
renormalization group calculation to determine exactly (within the limits of observation) 
the critical exponents which define the critical point in all dimension of phase transition 
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for the self organizing state of the system. Dynamic renormalization group calculation is a 
process which describes the relationship between levels in a multilevel system, such as 
atoms, molecules, cells, organs, bodies. A reference which explains how to calculate the 
critical point is the paper of Hwa and Carter, "Dissipative Transport in Open Systems: 
5 An Investigation of Self-Organized Criticality," Physical Review Letters , Vol. 62, No. 16 
(17 April 1989) pp. 1813-1816, the content of which is incorporated by reference. 
However, the subject matter is not considered an element of this invention. 

The simulation technique according to the invention will yield a 
confirmation of the prediction of the inherent critical point of each closed node whereby a 
10 stable system will result. Alternatively it will yield sufficient information to evaluate the 
deviation from the stationary state in the subject system. For example, in a physiological 
system, deviations from steady state are deviations from homeostasis which correspond to 
illnesses. 

The invention has been explained with respect to specific embodiments. 
15 Other embodiments will be apparent to those of ordinary skill in the art. A technique has 
been disclosed for simulating complex systems in terms of simple sources, simple sinks 
and simple nodes with critical points scalable across dimensions. Thus, this invention is 
not limited, except as indicated by the appended claims. 



WHAT IS CLAIMED IS: 
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1 1. In a computer system, a method for simulating a dynamic system 

2 with a plurality of interacting nodes of interest in a network of said nodes of interest, said 

3 method comprising: 

4 providing said nodes of interest in said computer system, each node of 

5 interest having at least one input, at least one output paired with said at least one input, at 

6 least one transformation of inputs, at least one transformation of outputs, a measurable 

7 ratio of input transformation rate to output transformation rate of an input/output pair, at 

8 least a first activated state in the node corresponding to an excess measurable ratio of 

9 input to output, at least a second activated state in the node corresponding to a deficit 

1 0 measurable ratio of input to output, and transient storage of a product of the input; and 

1 1 for each node of interest, defining a balanced state between said first 

12 activated state and said second activated state, said balanced state corresponding to a zero 

13 error between said measurable ratio and a preestablished balanced ratio, said 

14 preestablished balanced ratio corresponding to a mathematical critical point in 

15 thermodynamic energy. 

1 2. In the method according to claim 1 further including the steps of: 

2 for each said node of interest, sensing for non-zero error between said 

3 measurable ratio and said preestablished balanced ratio; and 

4 using said non-zero error as a control signal to mediate at least one of said 

5 inputs, said outputs and an external process. 

1 3. The method according to claim 2 wherein said node is 

2 representative of a living organism and wherein said error signal provides input to a 

3 regulating element for regulation to a condition of homeostasis. 

1 4. The method according to claim 2 wherein said node is 

2 representative of a non-living system and wherein said error signal is at least an indication 

3 of imbalance in energy distribution. 

4 Pathways span multiple elements in a system across multiple dimensions. 
1 5. The method according to claim 2, further including: 
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2 establishing pathways between outputs of first selected nodes of interest to 

3 inputs of second selected nodes of interest. 

1 6. The method according to claim 2 further including depicting each 

2 said four dimensional model in five orthogonal dimensions of space, time and grayscale, 

3 said grayscale representing a mapping from a second temporal dimension. 

1 7. The method according to claim 6 further including providing 

2 feedback across said five orthogonal dimensions from said old four dimensional model to 

3 produce a new four dimensional model, said old fo ur dimensional model and said new 

4 four dimensional model together constituting a six dimensional model. 

1 8. The method according to claim 1 wherein said critical point is 

2 selected for maximum stability of said balanced state. 

1 9. The method according to claim 1 wherein said critical point is 

2 selected in response to sensing said outputs of said nodes. 

1 10. In a computer system, a modeling node for use in simulating a 

2 dynamic system in a network of said nodes, said node comprising: 

3 at least one input; 

4 at least one output paired with said at least one input; 

5 at least one transformation of inputs; 

6 at least one transformation of outputs; 

7 a measurable ratio of input transformation rate to output transformation 

8 rate of an input/output pair; 

9 at least a first activated state in the node corresponding to an excess 

1 0 measurable ratio of input to output; 

1 1 at least a second activated state in the node corresponding to a deficit 

12 measurable ratio of input to output; 

13 transient storage of a product of the input; and 

14 a balanced state between said first activated state and said second activated 

15 state, said balanced state corresponding to a zero error between said measurable ratio and 

16 a preestablished balanced ratio, said preestablished balanced ratio corresponding to a 

17 mathematical critical point in thermodynamic energy. 
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METHOD AND APPARATUS FOR DESCRIBING AND SIMULATING COMPLEX 

SYSTEMS 
ABSTRACT OF THE DISCLOSURE 
A method for description and simulation based on organizing data into 
maps of invariants, the invariants being points of energy balance in a system of interest 
which is either in a stationary state or in a transitory disturbed state. The method includes 
identifying invariants in the system of interest by identifying primary sources and sinks of 
energy, identifying secondary energy sources and sinks coupled to the primary sources 
and sinks, and coupling all such sources and sinks into a network of transformations 
organized around nodes of those sources and sinks corresponding to the invariants, each 
of the nodes being characterized by a locally defined principle of balanced self- 
organization in a system with both a conservation law and energy dissipation. Such a 
system becomes "organized" upon achievement of a critical rate of entropy flux into the 
environment. Associated with each invariant are response rates related to energy transfer 
rates into and out of the invariants. The invariants are mathematically similar to the 
critical point found in equilbrium systems that undergo second order phase transitions. 
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